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Abstract 


A numerical simulation of radiative, conductive, and 
convective heat transfer of the Martian dust~laden atmos- 
phere-soil system is presented with particular emphasis^ 
given to heating/cooling in regions of sharp variation in 
temperature or absorption and its resultant impact on out- 
going planetary spectral radiance, as measured by the 
Mariner 9 IRIS. Thermal coupling betv;een the ground and 
atmospheric subsystems is modeled by the total heat flux 
balance at the interface. In the simulation procedure, 
local thermpdynamic equilibrium (LTE) is assumed, and a 
combined strong -weak line transmission function permits 
short- and long-range exchanges of energy from the surface 
toward space. Direct absorption of insolation in the near- 
IR bands by both silicate dust and CO^ is incorporated. 

The thermal coupling appears as a boundary- initial value 
problem which has been solved through forward -time inte- 
gration in 15-minute steps. The input data base stems from 
various experiments flown on Mariner 9. Results of interest 
are: (1) atmospheric counter-radiation is increased by more 

than 50% in the presence of dust, (2) tropospheric convective 
mixing is confined to the lowest 5 km, (3) quasi-isothermal 
stratification results in the tipper part of the dust layer 
during the period of intense solar heating, (4) a highly 
effective radiation-convection layer exists at the top of the 
dust stratum with a diurnal temperature variation of 40 K. 
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Exne numerxsche Simulation des V/arme transports 
durch Strahlung, Warmeleitung und Konvektion in dem 
ataubgeladenen Atmosphare-Boden System des Planeten 
Mars wird vorgelegt, die insbesondere Ewarmung und 
Abkuhlung in den Bereichen scharfer Anderung der 
Temperatur Oder Absorption berucksichtigt mit dem 
sich daraus ergebenden Einfluss auf die aufwarts 
gerichtete, planetarische , spectrale Radianz, wie sie 
durch Mariner 9 IRIS gemessen wurde. Die thermische 
Koppelung zwischen dem Boden und dem atmospharischen 
System ist modellmassig gegeben durch die Bedingung 
der Kontinuitat des totalen Warmeflusses an der 
Zwischenf lache . Der Modellierungsvorgang sieht lokales 
thermodynamisches Gleichgevicht (LTE) vor; eine 
kombinierte Transmissionsfunktion fur starke und 
schwache Absorptionslinien erlaubt Energieaustausch uber 
kurze und weite Distanzen vom Boden bis zu grosser Hohen- 
Ausserdem ist direkte Absorption der Sonnenstrahlung 
im nahen Infrarot einbezogen fur Silikatstaub und 
Kohlendioxyd . Die thermische Koppelung stellt sich als 
Grenz-Anf angswert-Problem dar , das durch Vorwartsintegration 
in 15-Minuten Intervallen gelost worden ist. Die 
Eingabedaten stanunen von verschiedenen Experimenten , die 
mit Mariner 9 geflogen wurden. Ergebnisse von Interesse 
sind: (1) atmospharische^Gegenstrahlung ninunt urn mehr als 

50% zu, wenn Staub gegenwartig ist; (2) tropospharisch- 
konvektive Mischungsvorgange sind auf die untersten 5 km 
beschrankt; (3) quasi-isothermische Schichtung ergibt sich 
im oberen Teil der Staubschicht w^hrend der Zeitspanne 
intensiver solarer Erw^mung; (4) eine Schicht hochwirksamer 
Strahlungsdivergenz und Konvektion existiert an der 
Obergrenze der Staubschicht, die eine ganztagige Temperatur- 
schwankung von etwa 40°K aufweist. 



1. Introduction 


The global extent of sand-like dry soil material on Mars 
permits the surface temperature field to respond rapidly to 
changes of largely-unattenuated insolation , with the mid- 
latitudinal diurnal variation being ~70K. In addition, the. 
predominance of CO^ and relative tenuity of the .Martian 
atmosphere result in a large radiative response capability; 
therefore, a significant diurnal thermal boundary layer is 
induced even if convection were absent. Analysis of radiative 
relaxation times for Mars (Goody and Belton, 1967) indicates 
a characteristic depth of 1.6 km for the diurnal temperature 
wave. The latter study also shows that molecular conduction 
may be an important transfer mechanism for. thermal per- 
turbations of spatial scale < 0 (lO.m ) . 


The purpose of this note is to present some first findings 
of a numerical modeling developed to (1) investigate the 
thermal coupling of the Martian atmospheric and soil sub- 
systems , and the resultant diurnal temperature variation in 
the lowest scale- heights of the atmosphere and ( 2 ) develop 
some criteria for the interpretation, under the condition of 


atmospheric dust loading. 


of i-iarinsr 9 IRIS inverted data. 


The design of the numerical model departs in several 
essential ways from that employed by Gierasch and Goody (1967, 
19B8). In the latter model, a "traDotoidal rule" was used to 
approximate the exchange integral for planetary radiative 
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heating, correctly accounting for the boundary terms 
while representing the contribution from various atmos- 
pheric source layers with emission terms evaluated at 
two discrete points which are sometimes far apart. This 
approximation emphasizes the long-range exchange but 
eliminates the effects of short-path, strong-line radiative 
transfers. On the other hand, a strong-line approximation 
to the transmission function was used throughout. A bulk 
transmission for the 15jt CO 2 band was utilized, without 
temperature dependence in either the integrated line in- 
tensity or Lorentz half-width. As Gierasch and Goody have 
pointed out, the resulting model is not effective for 
phenomena involving small length scales., such as radiative 
heating/cooling near a region of sharp variation in tempera- 
ture or absorption. 

It is precisely the latter phenomena which are of greatest 
interest in the present study. VJe are especially concerned 
with (1) the diurnal history of thermal structure near the 
soil-atmospheric interface and the extremities- of a suspended 
dust layer, and with (2) the impact of relatively small-scale 
temperature variation on the outgoing Dlanetary SDectral 
radiance such as would be measured by the Mariner 9 Infrared 
Interferometer Spectrometer (IRIS-). The model described in 
the present paper is most effective for study of variations 
in heat exchange and thermal structure in the f ree-atmosphcric 
layers above the daytime convective boundarv laver, and in 
the near- surface regions V7hen the boundary laver is 
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stably-stratified. Expansion of the model to account for 
a time-dependent free-convective boundary layer will be de- 
scribed in a forthcoming article. In the present paper, we 
report these findings which are expected to be little affected 
by the inclusion of the convective boundary layer. 

2. Simulation Procedure 

The atmospheric model adopted is that of a 100% CO^ atmos- 
phere in local thermodynamic equilibrium (LTE) , with olane- 
parallel stratification and negligible Doppler broadening 
(Pallmann, 1968; Pallmann and Dannevik, 1972). The lower 
solid surface is assumed horizontal and of silicate composition. 
The effect on heat exchange of dust suspended in the lowest 
scale heights is incorporated by considering the increase of 
effective radiative ' path length induced by multiple scattering 

and the additional direct absorption by the dust particulates 
themselves. 


From dynamical scaling considerations and radiative re- 
laxation time estimates, Gierasch and Goody (1968) have in- 
ferred that in the first approximation, temperature change 
due to advection by the mean wind field may be neglected com- 


pared with that due to radiative, molecular-conductive, and 
turbulent-convective transport. Thus, the atmospheric response 
to surface temperature variation may be investigated without 
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Details of the coupling beti-zeen Martian soil and atmospheric 
temperature variation depend crucially on the structure of 
the CO^ absorption bands. Since most of the 15/^ -band 
absorbs strongly even under the reduced surface pressure 
conditions of Mars 6 mrb) , a portion of the coupling will 
occur through the near-surface layers. Strong-line radiative 
transfer will approximate a diffusion process in these layers , 
augmented by the conductive-convective heat transports induced 
by large temperature gradients. On the other hand, spectral 
"gappiness" and weak-line transmission, each enhanced by 
atmospheric tenuity, will permit some long-range exchange of 
energy from the surface to several scale heights and toward 
space (Drayson, 1972). This longer-range exchange approxi- 
mates a Newtonian cooling (Kuo,. 1963 ; Schlichting, 1960). 

Several factors serve to complicate the integrated radi- 
ative— conductive-convective exchange process. Among these are 
direct absorption of insolation in the^near-IR CO^ bands, the 
presence of spectral lines of intermediate intensity, sub- 
stantial variations in gas-kinetic properties along the radi- 
ative path length, and the presence of dust-aerosol in the 
troposphere. One of the effects of the latter is to increase 
gaseous absorption, since multiple scattering results in a 
longer effective radiative path length. In addition, the dust 
□articulates may directly absorb both solar and planetary 
radiation. 
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At the soil-atmosphere interface, a discontinuity in 
physico-optical properties occurs, although on a physical 
basis we may still expect the resultant net flux of heat 
to remain continuous across the interface* The net heat flux 
includes contributions from short- and long-wave radiation 
and molecular-conductive fluxes in the soil and atmosphere. 

An additional physical constraint is a tendency tov/ard a 
matching of temperature determined by molecular conduction 
in the soil sub- system and the radiative-conductive temper- 
ature in the atmospheric layers. This tendency is comple- 
mentary to the LTE condition. 

In devising the simulation procedure, we have incorpo- 
rated to some exrent each of the above mentioned physical 
processes into a computer model, except for turbulent- 
convective heat transport. To accommodate the latter, the 
computer-generated temperature profiles, have been submitted 
to a convective adjustment in the appropriate atmospheric 
region above the constant flux layer,. The depth of this . 
convective region has been determined on the basis of com- 
plete adiabatic mixin.g. 

Thermal coupling of the soil and atmosphere may be 
studied as a boundary- initial-value problem through a 
forward time integration of the applicable heating rate 
equation. In our context, this relation takes the form. 
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f oo ^ 

K^|yS^(0)exp(-J'pK^dsu'‘ ) + 

X5 X 

G^(t)exp(-b/pK^ds) + J ( C, t )^^exp(-b Tp K^ds )J 

+ JBv(C,t )^[exp(-bj”pic^ds)] dc 2B^(x,t)|dv, 


0<x<Xs, t>0 


(la) 


PsCsl?2(x,t) = ^(ks|Is) 


3t 


X3<X<X<^, t>0 


(lb) 


where : 


T,p,c,k temperature, density, specific heat, and 
thermal conductivity, .in the atmospheric 
(no subscript) and soil subsystem (sub- 
script "s") 


x,C 

Xs 

Xd 

arccos ( p) 


depth, measured from top of atmosphere 

soil-atmosphere Interface level 

lowest level considered in soil subsystem 

solar zenith angle at time t 

mass absorption coefficient at frequency v 


B y 


specific intensity of blackbody emission 


upwardly-directed specific intensity at 
the interface level due to thermal emis- 
sion or diffuse reflection of solar irradi 
ation 

spectral solar radiation at the actual 
Sun -Mars distance 


b 


diffusivlty factor. 



The Tunc t ions and depend on depth and time, 

while u and depend on time. 


Continuity of the total net flux of heat at the 
interface requires that 

Pn(xs-0,t) - Pn(xs+0,t) = 0. (2) 


In terms of the relevant component fluxes, this re- 
lation represents an Internal boundary condition for 
the level x^: 



aeT^(xs-0,t) + 


kl^ + kglZs 

3x ^3x 


0 




(3a) 


Here, o, e, and represent the Stefan-Boltzmann 
constant, IR-emisslvity, and surface albedo. Match- 
ing of temperature at the interface is expressed as 


T(xs-0,t) = T(Xs+0,t). 


(3b) 


In (1) and (3a), it is assumed that no phase changes 


of CO 2 occur. 



The upper boundary condition at an elevation of 50 km 
(x-0) is set as vanishing incoming thermal radiation and 
constant temperature. Finally, the lov/er boundary condition 
at is that of constat temperature* 

A quasi— random model using the band parameters of 
Prabhakara and Hogan (1965) has been chosen for the trans- 
mission functions appearing in (la) and (3a), and a weak- or 
strong-line limit is employed, based on the magnitude of 
radiative path length pertinent to a given spectral and 
geometric interval* Frequency integration is performed over 
62 intervals embracing the 1-6^ and 12-l8^ vibration-rotation 
bands. In those atmospheric strata assumed dust-laden, the 
gaseous transmission is reduced by a factor exp(- 3 (i&x)^ where 

is an "effective absorption" coefficient due to dust, and 
Ax is the thickness of the layer. The cumulative fraction 
of irradiated solar energy flux v<hich is contained in the 
modeled CO 2 near-IR bands constitutes about 9^ of the total 
Insolation, With = 0,1km and the total depth of the 
dust layer being 30 km, the model thus produces approximately 
an additional 9^~attenuation of the vertical solar radiation, 
v;hich Is converted into heat v/lthin the dust layer* Gierasch 
and Goody (1972) utilized a 10^-attenuation under grey ab- 
sorption. 



To accomodate the inhomogeneity of temperature, 
pressure, and density, the transmission function for an 
extended path is approximated by a product of transmission 
functions taken over a number of intervening sub— layers, 
within each of which the gas-kinetic properties do not vary 
appreciably. The dependence of line intensity and Lorentz 
half-width on temperature and effective pressure is included 
in the calculation of individual layer transmission functions 

The integrals over depth in (la) and (3a) must be eval- 
uated by numerical quadrature. In view of spectral gappiness 
in the absorption spectrum and inclusion of molecular thermal 
conduction in the heating-rate equation, v/e anticipate a 
relatively complex structure of the vertical temperature pro- 
file in the boundary layer. A 52-point quadrature formula 
has been adopted, with 7 points distributed in the lowest 
km, and generally a 1 km spacing in the region between 1 
and 50 km. This resolution restricts the relative variation 
of Planck function to 3^ in a given layer, under rep- 
resentative conditions. The same spacing is used to con- 
struct finite difference approximations to the spatial 
derivatives in (la), while a 2-cm increment is used in the 
soil sub-system. 



Prom the estimates of Goody and Belton (I 967 ), the 
characteristic time for relaxation of thermal perturbations 
S order of several meters is in the range 

of lo3 sec for the Martian troposphere. On these physical 
grounds, we may infer that a computational time Increment 
of similar size is required to adequately resolve in time 
the propagation of temperature waves in the near— surface 
region. A 15-minute time step was used for all calculations 
discussed in this note. 

The computational procedure begins with calculation of 
transmission functions, solar spectral irradiance, and 
atmospheric spectral emission corresponding to the initial 
time and temperature distribution* These quantities are 
then combined to form the total (radiative + conductive) 
heating/ cooling function represented on the right-hand- 
side of (la). Flux components in (3a) are then computed, 
and this relation and (lb) are solved "for the soil tem- 
perature field for the end of the first timestep. Then, 

(la) is solved for the new atmospheric temperature field, 
using (3b) as the lower boundary condition. The new soil- 
atmosphere temperature profile is then used as input for 
the next timestep, and new solar spectral irradiances are 
calculated for the advanced local time. In this manner 



the procedure continues through the total interval of 
time desired for the simulation. 

Numerical values for physico-optical properties of the 
simulated soil material were derived through an Independent 
parametric study, in which the JPL (1968) thermal inertia 
values were adjusted until calculated diurnal surface 
temperature variation matched that suggested by recent 
Mariner radiometric scans. Raw spectral solar radiation 
data are those reported by Robinson ( 1965 ), adjusted to the 
Martian season and latitude relevant to a given simulation. 

The initial atmospheric sounding data. are derived from 
preliminary Mariner 9 IRIS inversion profiles kindly 
furnished by Drs. Hanel and Pearl of the Goddard Space Plight 
Center, and correspond to 1900 local Martian time (LMT), 
latitude 38 ®S, and solar declination r 23 ° (Southern hemi- 
sphere summer). ^ 

3 * Some results and preliminary interpretation 

Numerical output of the model for each timestep basi- 
cally consists of the thermal structure betv;een the levels 
of 60 cm soil depth and 50 km atmospheric height, radiative 
fluxes of solar, atmospheric, and surface origin in each of 
62 wave number intervals betv/een 1*29 and l8,02^, and the 



radiative and molecular-conductive heating/cooling rates 
at each computational level. Results presented here are 
based on two simulations, each representing an interval 
of 24 Martian hours. Jhe first case, designated "dust- 
free , utilizes as initial data the Mariner 9 IRIS in- 
version profile (Hanel et al., 1972) mentioned previously, 
with a starting time of 1900 LHT. The second case uses 
an identical initial sounding and starting tine, but 
includes a simulated ground-based dust layer of 30 )cm 
thickness; this case will be termed "dust-laden". 

Since over the night-time hours the radiative con- 
ductive model is thermally dissipative, minor incon- 
sistencies between the arbitrary subsoil and surface tempera 
ture and the atmospheric sounding tend to dissappear with 
time, so that by sunrise, the entire temperature profile 
is thermally "initialized". 

Fig. 1 depicts modeled atmospheric soundings for 0500 
and 1600 LMT for each simulation. With regard to the dust- 
free case, three basic regimes are evident. The lowest 
regime, from the surface to about 5 km, includes those 
atmospheric strata dominated by strong-line, short-path, 
diffusion-like radiative and convective exchange associated 
with soil surface temperature variation. The regime between 



5 and 31 km is "in communication” with surface tempera- 
ture variation primarily through weak-line exchange, but 
also undergoes direct solar heating and weak-line losses 
to outer space. The regime above 31 km is dominated by the 
upper boundary conditions , i .e . , direct solar heating and 
emission to outer space, and is affected little by soil 
surface temperature variation, which lags solar forcing by 
roughly 4 hours. 

Soundings from the dust-laden case exhibit essential 
differences from those of the dust— free simulation primar- 
ily in the middle regime (5-31 km; placement of the top of 
the simulated dust layer at about the same level as the 
upper bound of the middle regime for the dust-free case is 
coincidental). Over the daylight hours, this region 
approaches an isothermal stratirication , in agreement with 
recent Mariner 9 findings showing a tendency towards this 
^f^ucture in the dust— laden mission .^phases . The top of the 
dust layer acts as an "effective radiation surface", similar 
to its radiative response to the solid soil surface. Diurnal 
temperature variation at the top of the dust layer is nearly 
50% of the amplitude of surface temperature variation. Thus, 
the upper regime (above 31 km) reacts to this effective 

^<3.ble in the same v/ay that ihe lowest regime acts 
to surface temperature variation in the dust-free case. In 
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particular, the radiation table induces a statically un- 
stable stratification at the top of the dust layer in the 
afternoon hours. The lapse rate of the layer was subse- 
quently adjusted to the adiabatic value to reflect convective 

• • 

mixing. 

The Mariner 9 IRIS profile was derived by inverting 
the IRIS radiance measurements obtained during a dust-laden 
mission phase, to generate a temperature profile assuming 
a dust-free atmosphere. This initial temperature profile 
approximates a smoothed version of the modeled 1600 LMT 
output for the dust-laden simulation, with decreasing fidelity 
towards the surface. Further study is required to determine 
to what extent this result is coincidental or significant. 

As indicated on the right-hand-side of Fig. 1, the dust- 
laden profile at 0500 LMT is warmer than the dust-free in 
the lowest 800 m, an indication of the insulating effect of 
the dust blanket keeping surface temperature slightly higher 
in the nocturnal hours. However, in the bulk of the dust 
layer, a cooling has occurred by 0500 LMT which is sub- 
stantially larger than in the dust-free case. This is thought 
to be attributable to the adopted ’’grey-body" emission of 
the dust particles, V7hich results in more photons generated 
in spectral intervals outside the CO^ absorption bands. 
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Atmospheric and subsoil temperature profiles in the 

layers near the interface are shown in Fi^. 2, for the dust- 

free case. The early-morning atmospheric profiles exhibit 

diffusionlike heat transport characteristics, due to the dual 

action of strong-line railiative and molecular-conductive 

transfer. An "effective radiative-conductive diffusivity" 

2 

may be estimated by L /f v;here L andT are the characteristic 

length and time scales of the Diffusion" process. This 

5 2 —1 

diffusivity is found to be ^ 2 x 10 cm sec , on the basis 

of data represented in Fig. 2. This compares well with 

5 2 

Gierasch and Goody *s (1968) scaling estimate of 10 cm 
sec and is similar in order-of-magnitude to. a free- 
convective eddy diffusivity under mildly-unstable conditions. 
It also tends to corroborate Goody and Belton's (1967) 
estimate that the time for relaxation of thermal perturba- 
tions by radiation and mild turbulent convection should be 
comparable for a perturbation length scale of hundreds of 
meters. 

Each of the various heat flux components acting at the 
soil-atmosphere interface are shown as a function of the 
time of the day in Fig. 3. The limited computational resolu- 
tion in the atmospheric surface boundary layer cannot produce 
temperature gradients steep enough to model realistically 
the removal of heat from the interface by atmosoheric con- 
duction. fiowever, the MET curve in Fig. 3 indicates the 
energy surplus or deficit available 


to change the interface 
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temperature. Short-term oscillations of decreasing ampli- 
tude evident in the first two hours of the simulation 
represent the dissipative smoothing of inconsistencies in 
the arbitrary initial temperature distribution in the inter- 

I* 

face and subsoil layers. 

The fact that the model does not reproduce the initial 
flux values after 24 hours can be traced to at least two 
possible causes. One cause related to limited computational 
resolution has been discussed in the previous paragraph. A 
second possible cause is that the initial atmospheric tempera- 
ture profile may not be representative of the local time and 
season relevant to the IRIS measurement on which the inverted 
sounding is based. The preliminary Mariner 9 sounding which 
we utilize was based on an inversion method which assumes a 
dust-free atmosphere. 

Molecular-conductive heat flux in the soil material is a 

significant fraction of the net flux available for temperature 

* 

changes, at all times during the 24 hr simulated period. The 
conductive flux reaches a maximum about 1.5 hours before local 
noon, in agreement vjith measurements under comparable 
terrestrial physical settings (Sellers, 1965). 

Atmospheric counter-radiation varies little over the 
diurnal cycle. Since the concentration distribution of radia- 
t ively-active gases remains constant in the model, ary 
temporal variation in counter- radiation can be traced directly 
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to variations in the temperature distribution, for the dust- 
free case. The addition of a simulated dust layer results 
in an increased counter-radiation for all times of the day , 
as evidenced in Fig. 3. This phenomenon accounts for a 
higher minimum morning temperature of 19 8 K for the dust-laden 
case, as compared with 19 3K for the dust-free simulation. 

Height-time cross-sections of the temperature field be- 
tween 4 and 18 km of altitude are shown in Fig. 4a, b. Be- 
cause the pattern between the surface and 4 km is nearly un- 
changed by the presence of dust, it is not shown here. The 
region between 4 and 18 km is drastically modified by the 
dust layer during the daylight hours. This regime responds 
to direct solar heating and to a lesser extent, to interface- 
temperature variation. The dust-laden layers above 7 km 
develop a nearly isothermal stratification in the afternoon 
hours, in agreement with recent Mariner 9 observations. 
Detailed interpretation of additional features in Fig. 4 will 
be presented in the near future. 
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FIGURE LEGENDS 

Pig. la. Modeled vertical temperature profiles at 0500 
and 1600 local Martian tine (UdT), for the dust-free (solid 
) ^nd dust— laden (dashed line) conditions. Also shown 
is the initial (1900 LMT) profile obtained from inversion 
(Hanel, et al., 1972) of Mariner 9 IRIS measurements by 
assuming a dust— free atmosphere. 

Fig. lb. Profiles of the difference AT = in 

C Q 

temperature between the dust-free (subscript : c) and dust- 
laden (subscript ; d) simulations, at 0500 and 16OO liiT. 

All results based on follov;ing simulation parameters: 
latitude = -38^, solar declination 6 = - 23 ° (Southern 
hemisphere summer), surface albedo - 0.25, surface 
infrared emissivity e = 1,0, soil thermal conductivity 
ks = 2.3 X 1 q 5 erg cm ^ sec”^ K“^ (low-porosity silica),^ 
and C02-gas thermal conductivity k *= 1.33 x 10^ erg cm”^ 
sec“^ K*“^- 

Fig, 2 . Vertical temperature distributions in the sub- 
soil and lower atmospheric layers for selected hours of 
the Martian day, from the dust-free simulation. For 
environmental parameters used, see Fig. 1 , 
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Fig. 3. Diurnal history of the various component heat 
fluxes at the soil-atmosphere interface. EMIS represents 
the flux of long-wave radiation emitted by the surface, 

SCON the molecular-conductive heat flux in the soil 
material, SOL the absorbed insolation, ACR the absorbed 
flux of atmospheric counter-radiation, and NET the alge- 
braic sum of all component heat fluxes. Solid lines depict 
results for dust -free conditions, while the dashed lines 
(shown for EMIS, SOL, and ACR only) relate to the 
simulated dust-laden conditions. Influence of C02-gas 
thermal conduction is of lesser significance at the inter- 
face, and curve is not shown. 

Fig. 4a, b. Height-time cross-sections of the modeled 
temperature field (deg. K) between 4 and 18 km for (a) 
dust-free, and (b) dust -laden conditions. Approximate 
extremal lines are dashed. 
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